Condensation of magnons and spinons in a frustrated ladder 
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Motivated by the ever-increasing experimental eflort devoted to the properties of frustrated quan- 
tum magnets in a magnetic field, we present a careful and detailed theoretical analysis of a one- 
dimensional version of this problem, a frustrated ladder with a magnetization plateau at m = 1/2. 
We show that even for purely isotropic Heisenberg interactions, the magnetization curve exhibits a 
rather complex behavior that can be fully accounted for in terms of simple elementary excitations. 
The introduction of anisotropic interactions (e.g., Dzyaloshinskii-Moriya interactions) modifies sig- 
nificantly the picture and reveals an essential difference between integer and fractional plateaux. In 
particular, anisotropic interactions generically open a gap in the region between the plateaux, but 
we show that this gap closes upon entering fractional plateaux. All of these conclusions, based on 
analytical arguments, are supported by extensive Density Matrix Renormalization Group calcula- 
tions. 



I. INTRODUCTION 

Quantum phase transitions- occur when the variation 
of an external parameter (pressure, chemical composition 
etc.) produces a singular change of a system's ground 
state. If the transition is continuous, the properties of the 
system in the vicinity of the quantum critical point are 
dominated by universal features, which can be described 
in field-theoretic languagei Quantum phase transitions 
and related critical behavior are observed in a variety 
of physical systems, from heavy-fermion compoundaS, to 
quantum magnets^ and cold atomic gases^ 

In antiferromagnets the external parameter is typically 
an applied magnetic field H — {0,0, H). If the longitu- 
dinal component of the total spin, S"^, is a good quan- 
tum number, the magnetic moment of the ground state 
M = gusiS^) as a function of H may exhibit several 
plateaux, on which dM/dH = 0, separated by regions 
of continuously varying magnetization. The ends of a 
plateau are quantum critical points separating an incom- 
pressible ground state with an energy gap from a com- 
pressible state with gapless excitations. Affleck'^ noted 
in the context of spin 1 chains that such a phase transi- 
tion is similar to the condensation in a system of interact- 
ing bosons, a point of view reemphasized and extended 
in the context of coupled spin 1/2 laddersi^ The mag- 
netic field and magnetic moment play the roles of the 



chemical potential and particle number. The condens- 
ing objects are magnons, the quasiparticles carrying spin 
AS^ = ±1. The field theory describing the bosons near 
the Bose-condensation point^ is equally applicable to the 
end points of a magnetic plateau. 

Not every magnetization plateau ends in a simple con- 
densation of magnons. If the plateau state breaks a sym- 
metry of the lattice while the gapless state does not, the 
transition must restore the broken lattice symmetry. The 
universal properties of such a transition are expected to 
be different. In one spatial dimension, the picture based 
on the magnon condensation is generally applicable to 
"integer" magnetization plateaux defined as follows:^ the 
spin per unit cell differs from the maximal value by an 
integer. A fractional magnetization plateau may end in 
a phase transition belonging to a different universality 
class. Examples of such behavior were recently discussed 
by a number of authors j^iiSiii 

In this paper, we present a theoretical study of quan- 
tum phase transitions in a one-dimensional model anti- 
ferromagnet exhibiting both integer and fractional mag- 
netization plateaux. We employ numerical methods to 
observe critical behavior and compare the results to pre- 
dictions of the appropriate field theories. The outcome 
of our work stresses the importance of anisotropic in- 
teractions in the vicinity of quantum critical points, a 
point raised previously by previous authors^^ii^ Even a 
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FIG. 1: (Color online) The frustrated ladder. The exchange 
couplings are set to 1 on the vertical rungs, Ji — 0.55 on 
horizontal legs, and J2 = 0.7 on the diagonals. 

weak anisotropy makes a significant impact on the quan- 
tum phase transitions in question and, furthermore, the 
effects vary substantially between different universaUty 
classes. 

The paper is organized as follows. The model sys- 
tem and details of the numerical method are described 
in Sec. ^ In Sec. IIIII we discuss the ground states and 
phase transitions in the model with isotropic interactions. 
The influence of anisotropy is the subject of Sec. lIVI and 
a summary of the results is given in Sec. 



II. DESCRIPTION OF THE MODEL 

A. Hamiltonian and ground states 

Our model system is the frustrated ladder with spin- 
1/2 spins and with Heisenberg exchange!^ The largest 
exchange couplings, of strength 1, are on the rungs, with 
somewhat weaker couplings Ji along the legs and J2 
along plaquette diagonals, see Fig. The Hamiltonian 
is 

L 

'H = ^[S„4 • S„,2 - H{S^ ^ + 3^ 2)] 

n=l 

L-1 

) (1) 

n=l 
L-1 

+ J2 E](S„,1 • S„+i^2 + S„^2 • S„+i^i). 

n=l 

A ladder without diagonal links (J2 — 0) was examined 
in Ref. The point Ji = J2 is special: the spin of 
each rung is a conserved quantity and the ground state 
is known exactly. Conditions for the existence of a m = 
1/2 plateau in the general case were given in Ref. IT^ 

In this work, we study the case where the inter- 
rung couplings Ji and J2 are nonvanishing and similar 
in strength. Proximity to the exactly solvable model 
Ji = J2 suggests a partition of the Hamiltonian into 
an exactly solvable part plus a small perturbation. To 
that end, it is convenient to introduce the operators of 
the total spin of a rung S„ — S„^i + S„^2 and the spin 
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FIG. 2: (Color online) (a) A sketch of the ground state at low 
magnetic fields: all of the rungs are in the S = state, (b) 
Elementary excitations out of this state are magnons carrying 
spin AS'' = -1-1. 



difference D„ — S„.i — S„^2- The exchange part of the 
Hamiltonian is then a sum of three terms: 

L 

Ho - 5](S„)V2-i?^^, 

L-1 

'H\ = J E^ S„ • S„+i, (2) 

n=l 

"^2 = — ^Dn-Dn+i, 

n=l 

where J = (Ji + J2)/2 and 5, J — Ji — ,12- The algebra of 
D operators is discussed in Appendix 1X1 

The physics of the model is rather simple when the 
energy scales are well separated, 

(5J<J<1. (3) 

The dominant term Tip favors a spin singlet on every 
rung at low fields H < 1 as depicted in Fig.[5fa). In high 
fields, H > \, the state of lowest energy is the S"^ = -f 1 
component of the triplet. (The other two components of 
the triplet are high-energy states in any magnetic field.) 
The all-singlet and all-triplet states are the origins of 
the two integer magnetization plateaux with the spin per 
rung m — / L — and 1, respectively. 

The next term Tii represents a repulsion between 
triplets on neighboring rungs. This repulsion does not al- 
low the triplets to condense on all rungs at once and intro- 
duces an intermediate, fractional magnetization plateau 
with a triplet on every other rung and spin per rung m = 
1/2. Accordingly, there are two ground states breaking 
the translational symmetry, illustrated in Figs. OJa) and 
(b), which we refer to as the Neel states. The fractional 
plateau exists in the range of fields 1 < H < 1-1-2 J. 

Finally, the smallest term Ti.2 endows the triplets with 
mobility: unlike the total spin of the rung , the spin 
difference D„ does not commute with the dominant term 
Ho- The triplets acquire a kinetic energy of the order SJ. 
As a result, the fractional and integer plateaux become 
separated by gapless phases with a constantly varying 
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(d) 

FIG. 3: (Color online) (a) and (b) depict the two ground 
states of the fractional plateau with translational order. 
Elementary excitations are domain walls carrying = —1/2 
at (c) the low-field end of the plateau and — +1/2 at (d) 
the high end. 



magnetization and with spin correlations that decay as a 
power of the distance between sites. 

The resulting Hamiltonian ||2Jl has an axial 0(2) sym- 
metry. However, in real materials that symmetry is 
often violated by small additional interactions induced 
by the relativistic spin-orbital coupling. In many cases 
such interactions act like a staggered transverse magnetic 
fieldii2ii^ We consider the influence of such additional 
terms in Sec. IIVI 



B. Elementary excitations 

At the end of the m — plateau, the elementary ex- 
citations are individual triplets carrying spin AS^ = +1 
[Fig-EIb)]. Similarly, low-energy excitations near the end 
of the m = +1 plateau are isolated singlets in the back- 
ground of S"^ = -f 1 triplets. They carry spin AS^ = — 1. 
We refer to both of these excitations as magnons. 

The low-energy excitations near the ends of the frac- 
tional plateau are domain walls with spin AS^ = —1/2 at 
the lower-field end and AS^ = +1/2 at the high-field end 
[FigS-EIc) and (d)]. They will be referred to as spinous. 

The quantum phase transitions between the plateaux 
and gapless phases are triggered by the condensation of 
these elementary excitations. 



C. Mapping onto an XXZ spin chain, hard-core 
bosons, and fermions 

In the regime described by Eq. each rung can be 
found either in the singlet or in the = -1-1 triplet state. 
Effectively, we can treat this as a system with only two 
states per site, using perturbation theory in SJ to define 
a Hamiltonian that acts on this reduced Hilbert space. 

By identifying the singlet and the 5^ = 1 triplet states 
with the spin-up and spin-down states, respectively, we 
map the ladder onto a spin-1/2 XXZ antiferromagnetic 
chain with an easy-axis anisotropy^ Defining the usual 
spin-1/2 matrices s^, and s^ acting on the reduced 
Hilbert space, one has 

L-l 
n=l 

L 

- i/cdgc {sl + sl)/2- H^^, J2 < ' (4) 

where 

j^^dJ + 0{dJ^), j, = J + 0{6J^), 
H^^,^H~1~J + 0{SJ^). (5) 

The additional magnetic field -ffodgc = J at the chain 
ends breaks the symmetry between the two Neel states of 
this effective spin chain. The edge field plays a role in the 
formation of the ground state at the fractional plateau, 
where singlets and — +1 triplets have comparable 
energies. 

The XXZ chain is gapped in the antiferromagnetic 
regime jz/\jx\ > 1 and i?xxz = 0, with zero magneti- 
zation in the ground state»i^ The presence of the energy 
gap means that the magnetization remains exactly zero in 

a finite range of fields — i/min < ^^xxz < -ffmin- At ±Hynin 

the energy gap vanishes and the spinous condense. (The 
corresponding fields in the ladder corresponds to the edge 
of the M = 1/2 plateau and are called Hcs and Hc2, re- 
spectively). As |i?xxz| is increased further, the ground 
state becomes a sea of spinous with a growing magneti- 
zation. The system in this regime is a Luttinger liquid 
with continously varying critical exponents.^® Finally, at 
i/xxz = ±-f^max the magnetization of the XXZ chain sat- 
urates (in the ladder, this corresponds to the saturation 
field Hc4 and to the field Hd at the end of the m=0 
plateau). The critical fields of the XXZ chain are known 
exactly iii 

(—1)" 

Hma.x = jz + \jx\, i?min = \jx\sinhg , 

'^-^ cosh no 

n— — oo 

(6) 

where coshg = jz/\jx\- 

The system can also be viewed as a one-dimensional 
hard-core Bose gas: the singlet state of a rung is mapped 
onto an empty site, and & = +1 triplet becomes an 
occupied site. The bosons have hopping amplitude 5J/2, 



strong nearest-neighbor repulsion J, and chemical poten- 
tial -ffxxz- 

In one dimension, the hard-core constraint can be re- 
moved by replacing the bosons with spinless fcrmions. 
This representation is particularly convenient when the 
Bose system is nearly empty (m <C 1) or nearly filled 
(1 — m ^ 1). In these limits, the short-range re- 
pulsion between the fermions is largely suppressed by 
the Pauli principle. Near the fractional plateau, when 
\m — 1/2| <g; 1, it is convenient to represent the domain 
walls as spinless fermions i-iS. 



D. Numerical work 

When the energy scales do not form the hierarchy of 
Eq. 0, one must resort to numerical methods. Nonethe- 
less, the general picture painted above remains largely 
intact and, furthermore, the critical behavior near the 
quantum phase transitions is expected to be univer- 
sal. To verify this, we have numerically determined the 
ground state, its magnetization, and the energy gap in a 
ladder with coupling constants Ji = 0.55 along the legs 
and J2 = 0.7 along the plaquette diagonals. The ground 
and first excited states were determined by the DMRG 
mcthodi^ We have worked with ladders with up to 200 
rungs, and have used the so-called "finite algorithm" ver- 
sion of the DMRG method. The use of open boundary 
conditions allows us to study oscillations of magnetiza- 
tion (5^) induced by the presence of the ends. 

For the model with isotropic interactions, we have car- 
ried out two sweeps and have retained m — 600 states 
in the system block. The typical weight of discarded 
density- matrix eigenvalues is of order 10^^^. We have 
performed a few calculations with 6 sweeps and m ~ 1000 
states; the energy difference was of the order of 10^^. We 
have calculated the ground state in each Sz sector in zero 
magnetic field for sizes up to iV = 400 and by comparing 
energies in a field have deduced the global ground state 
as a function of magnetic field. 

For the model with a staggered field, Sz is no longer a 
good quantum number, and we need to make a calcula- 
tion for each value of the magnetic field H . Fortunately, 
due to the opening of the gap, far fewer states are needed: 
carrying out two sweeps and keeping only m = 200 states 
for up to = 200 sites and m — 400 for larger systems, 
the typical discarded weight was of order 10~^^. 

The other subtle issue when performing DMRG sim- 
ulations is the choice of boundary conditions. Since one 
cannot access all excited states with the DMRG, choos- 
ing the appropriate boundary conditions can be crucial 
to obtaining relevant information about the excitation 
gap. 

In the case of an isotropic chain, for which the gap is 
not an issue, we have worked with open boundary con- 
ditions. As a result, in the m = 1/2 plateau the ladder 
with an even number of rungs L has a single spinon in 
the ground state (see section III.B.2 for details). The 
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FIG. 4; (Color online) Magnetization per rung m = {S^)/L. 



small magnetization step in the middle of the fractional 
plateau occurs when the spinon changes its spin from 
AS^ = -1/2 to +1/2. 

For systems with a staggered field, however, we have 
used asymmetric boundary conditions, imposing different 
values of the rung couplings at the edge: J = 1.8 for the 
first rung and J — 0.2 for the last one. These boundary 
conditions are necessary for the determination of the en- 
ergy gap in the formerly gapless regime between Hd and 
Hc2 [Fig-E]- With open boundary conditions, the system 
has two nearly degenerate ground states. Altering the 
end rungs as described above lifts this near-degeneracy 
and pushes one of the ground states well above the first 
excited state. 



III. LADDER WITH ISOTROPIC 
INTERACTIONS 

In this section, we present a detailed analysis of the 
properties of the frustrated ladder without anisotropic 
interactions. The basic properties of this system have 
been described beforeii^*i^*i& In particular, the magne- 
tization curve is expected to have plateaux at m = 0, 
TO = 1/2, and to = 1, to increase monotonically between 
the plateaux, and to have square-root singularities at all 
plateau ends. In practice, the magnetization curve dif- 
fers from the naive expectations in a number of ways. 
For example, the square-root singularity in the magneti- 
zation near the fractional plateau is barely detectable; a 
jump in magnetization is seen in the gapless region be- 
tween the m = 1/2 and m = 1 plateaux; the magnon 
mass near Hd differs significantly from the prediction of 
perturbation theory. We discuss the main features of the 
magnetization curve below. 
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A. Critical points 

The results for L = 52 and 152 rungs are shown m 
Fig. 21 As expected, the magnetization curve has three 
plateaux at m = 0, 1/2, and 1. The m = plateau 
ends at Hd — 0.806, the m = 1/2 plateau lies between 
Hc2 = 1.345 and i/^a = 1.974, and a fully polarized state 
with m = 1 begins at Hd = 2.500. A finite jump in the 
magnetization Am = 0.099 is observed in a compressible 
regime at Hd = 2.254. We next discuss the critical fields. 

1. Magnon condensation at Hd 

It is straightforward to calculate the end point of the 
fully polarized state m = 1. When the energy scales 
separate well, i.e., in the regime described by Eq. JSJ, 
the excitations with the lowest energies are singlets with 
the dispersion 

eo^oik) = H -1-2J + SJcosk. (7) 

The singlets are gapped in fields H > Ho,o = 1 + 2 J + 
\6J\ = 2.4 for our choice of couplings. However, the 
numerics show that the condensation occurs at a some- 
what higher field. This discrepancy can be traced to a 
poor separation of energy scales: J = 0.625 is not that 
small. As a result, the first excitations to condense are 
the S'^ — components of the triplet, whose energy dis- 
persion is 

ei,o{k) ^ H -2J + 2Jcosk, (8) 

and the condensation field Hi^q = 4J = 2.5, in perfect 
agreement with the numerics. 

2. Spinon condensation at Hc2 and Hcs 

To lowest order in 6 J, the end points of the fractional 
plateau m — 1/2 can be obtained in a similar way. The 
spinous carrying — ±1/2 have energies 

e-i/2(fc) = {H -l)/2 + SJcos2k + 0{Sj'^), (9) 
e+i/2(fc) = [l- H)/2 + J + SJcos2k + 0{SJ^). 

These spinous condense at 

Hc2 ~ 1 + 2\6J\ = 1.3 and (10) 
Hc3 ~ 1 + 2J-2|(5J| = 1.95, 

respectively, which are not that far off from the values 
Hc2 = 1.345 and H^s = 1.974 obtained from the DMRG. 

To improve the lowest-order estimate, we have ex- 
panded the parameters of the XXZ chain ||SJ) and its 
critical field ffmin ® to 0{6J^). The physical origin 
of these corrections can be traced to quantum fluctua- 
tions of spins in the ground and excited states of the 
m = 1/2 plateau and its excitations. The term 7^2 in the 



Hamiltonian connects the = —1/2 spinon [two adja- 
cent singlets, see Fig.|3Ic)] to high-energy states, shifting 
its energy by an amount C 5,P < 0. The energy of the 
= -1-1/2 spinon [Fig. |3|[d)] is unaffected by quantum 
fluctuations at this order. 

The energy shift of the S"^ = —1/2 spinon can be taken 
into account by adding the following term to the Hamil- 
tonian of the XXZ chain: 

AHxxz = ^(l/2-<)(l/2-<+i)CW2 (11) 

= const -CSJ'Y. < + CSJ'Y. «+i ■ 

It can be seen that the added term affects both the Ising 
coupling and the effective field of the XXZ chain: 

j.,^dJ + 0{5J^), 3, = J + CSJ^ + 0{SJ^), 
H,,^r,^H-l~J + C6j'^ + 0{6J^). (12) 

Expansion of Eq. (jSJ in powers of ^ J yields the lower 
critical field of the XXZ chain 

Hmin = J- 2\6J\ + {C + 1/2 J)SJ^ + 0{SJ^). (13) 

The resulting condensation point of the = +1/2 
spinous 

i/c3 ~ 1 + 2J - 2|W| + (5JV2J = 1.968 (14) 

is now very close to the DMRG value of 1.974. 

The critical field of the 5^ = ^1/2 spinous is sensitive 
to the energy shift C5,P . The function C{J) is computed 
in Appendix m For J = 0.625 we obtain C ^ -2.026 
and thus 

i?c2 « 1 + 2\5J\ ~ {2C + 1/2 J)(5J2 = 1.374. (15) 

Comparing it to the DMRG value of 1.345, we see only 
a modest improvement over the first-order result HlOfl. 
The apparent reason for the slow convergence of the per- 
turbation series for Hc2 is the fairly small energy gap 
(A = 0.255) separating the 5*^ = —1/2 spinous from 
higher-energy states (see Appendix 0. 

3. Magnetization jump at Hd ~ 2.254 

Between Hc3 and Hd, the ground state switches from 
a mixture of S*^ = -1-1 triplets and singlets to one of — 
+1 and = triplets. This transition is accidental in 
the sense that it is not accompanied by any change in 
symmetry. It is therefore not surprising that the change 
is accompanied by a discontinuity in the magnetization. 

4-. Magnon condensation at Hd 

The energy of an isolated AS*^ = 1 excitation at the 
m = plateau is 

ei.iik) ^1- H + 5J cosk + 0{6J^). (16) 
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(a) 






FIG. 5: Corrections to the kinetic energy of a magnon at 
the order SJ^: (a)^(b) A pair of triplet states with net spin 
zero is created adjacent to an existing magnon. (b)-^(c) One 
of those triplet states annihilates with the original magnon, 
leaving a new magnon on a different lattice site. This process 
provides an additional channel for magnon motion and thus 
lowers the effective magnon mass. 



Thus, to first order in SJ, the magnon condensation is 
expected at TJ^i = 1 — \SJ\ = 0.85, not very far from the 
DMRG result 0.806. 

At 0{6J^) quantum fluctuations not only lower the 
energy of the magnon but also change its hopping am- 
plitude (Fig.l^J. The magnitude of the second-order cor- 
rection is rather large, again thanks to a small energy 
gap (A = 0.125) between the magnon and higher-energy 
states. See Appendix ICl for details. 



B. Magnetization patterns 
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FIG. 6: (Color online) Distribution of spin (5^) as a func- 
tion of rung position n in the ground state off of the m = 
plateau. The symbols are the numerical results at = 2, 
4, and 8. The curves are spin distributions for 2, 4, and 8 
magnons treated as free fermions carrying spin AS^ = -1-1. 
Calculations are on a ladder with L = 160 rungs. 
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Thanks to the presence of an energy gap, both the total 
spin and the average spin of an individual rung (S^) 
remain exactly zero in the low- field regime \H\ < Hd- 
As the field increases beyond Hd , both the local and the 
total spin begin to increase. In a finite ladder, the local 
magnetization (5',^^) is distributed in a nonuniform way, 
revealing interference patterns. 

1. m I: dilute magnons 

At low concentrations, the magnetization is carried by 
individual magnons which can be viewed as hard-core 
bosons or, more conveniently, as free fermions with spin 
AS^ = +1. (The nearest-neighbor repulsive potential 
acting between magnons is rendered irrelevant by the 
Fermi statistics.) Treating the magnons as ideal fermions 
one obtains a magnetization distribution 

fc=i 



FIG. 7: (Color online) The square of magnetization density 
m off the m — plateau. 

where, in the continuum approximation, ipkin) = 
\/2/L sin {-nkn/ L) is the wavefunction of a nonrclativis- 
tic fermion in a box of length L + I. This simple 
model agrees well with the magnetization distribution ob- 
tained numerically at small values of the magnetization 
m — S^/L (Fig. inj. Deviations already become notice- 
able when m reaches 1/20. 

The free-fermion approach predicts a fast initial 
growth of magnetization m = S^' / L: 

m - AS*" kp/TT = TT^i lAS-'l^/^ ^2m*\H - Hd\, (18) 

where the inverse mass l/m* w We have found 

numerically that indeed rises linearly with H (Fig.[7J). 
However, the slope dm? /dH is less than a third of the 
calculated value of 2tt^'^\5J\^^ . 

To track down the source of the discrepancy, we have 
computed the contribution of higher-order processes to 
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FIG. 8: (Color online) Distribution of magnetization {S'^ in 
a ladder with L = 160 rungs (a) Exactly at the magnetiza- 
tion plateau m = 1/2, or = I//2, the ladder contains one 
spinon. (b) At = L/2 — 1, three spinous are present. 

the magnon dispersion. The 0{5J'^) term turns out to 
be larger than the first-order result. This can be un- 
derstood on a qualitative level by considering a typical 
0{5J^) contribution to the kinetic energy of the magnon 
in which two additional magnons are created and de- 
stroyed (Fig. ISjl. The lowest-lying 3-magnon state has 
an energy of A = 2 — 3 J above that of a single magnon. 
For J = 0.625, this energy gap A — 0.125 is compara- 
ble to the small parameter 5 J — —0.15, meaning that 
the lowest-order perturbation theory in 5 J may not give 
reliable results. See Appendix lUl for details. 



2. \m — 1/2| ^ 1.- dilute spinous 

The numerically determined distribution of magneti- 
zation (S*^) exactly at m = 1/2 {S^ = 80 in a ladder 
with L — 160 rungs) is shown in Fig.lS^a). Contrary to 
our initial expectations, the plateau state is not a sim- 
ple Neel state with a constant staggered magnetization 
on top of a constant background 1/2-1- (— 1)"/2. The 
staggered magnetization is highly nonuniform. 

To understand this result, consider the m — 1/2 




(d) 

FIG. 9: Low-energy states of a ladder with an even number 
of rungs at magnetization density m = 1/2. 

plateau state in a finite ladder with an even number of 
rungs L. To ©(AJ^), its ground states are configura- 
tions with L/2 triplets with no two triplets next to each 
other. There are only two such states in a ladder with pe- 
riodic boundary conditions. In contrast, for open bound- 
ary conditions there are L/2 -I- 1 such configurations: two 
Neel states and L/2 — 1 states with a singlet-singlet do- 
main wall (Fig.O. The term 7^2 delocalizes the domain 
wall and is ultimately responsible for the strong modula- 
tion of the staggered magnetization. 

The mapping onto the XXZ chain provides an alter- 
native perspective. Ordinarily, the two Neel states of 
the XXZ chain are degenerate even in the presence of 
a uniform magnetic field. However, our system has 
an extra magnetic field of strength — J/2 at the ends, so 
that both end spins prefer the -1-1/2 state. In a chain of 
even length, this inevitably leads to the formation of a 
domain wall. The cost of the wall (-1- J/2) is exactly offset 
by the reduction of the energy of the end spin (—J/2); a 
negative delocalization energy of the domain wall (— 1(5 J|) 
lowers the energy of this state relative to the Neel state. 
Thus the effective XXZ chain always has spins -1-1/2 at 
the ends in this regime. Accordingly, the end rungs of 
the ladder are in the triplet state and thus there must be 
a single domain wall somewhere in the chain. 

The strong modulation of the staggered magnetization 
evident in Fig. |HIa) can be traced to the delocalization 
of the spinon. Rung n is in the state with spin S = 
1/2 — (— l)"/2if the domain wall is on its right, otherwise 
it has spin S = 1/2 -f (— 1)"/2. For a rung near the left 



edge, the domain wall is almost always on the right and 
vice versa. Thus we expect 



1 



(19) 



where the smooth envelope E{n) interpolates between +1 
and —1 and has a node in the middle. 

To evaluate the envelope E{n), we adopt the contin- 
uum approximation and treat the spinon as a nonrela- 
tivistic particle in a one-dimensional box. Doing so yields 

X / X I , / m9 , , 2n sin(27rn/i) 
E{n)= sgxv{x-n)\iPi{x)Y dx = \ — -+ ^ 

Jo L TT 

(20) 

where i/jiin) = ^2/Lsin {nn/L) is the ground-state 
wavefunction of a spinon. The result is shown in Fig.|SJa) 
as a dashed curve. While the agreement is already rather 
good, further improvements can be made, as explained 
below. 

Just to the left of the m = 1/2 plateau, the ladder con- 
tains a few spinons with spin A5'^ = — 1/2 each. Because 
the first and last rungs still remain in the triplet state, 
the ground state of a ladder with an even number of rungs 
L and spin 5'^ = L/2 — p contains r — 2p + 1 spinons. 
At low concentrations r/L the spinons can be treated as 
ideal fermions»iS The envelope of the staggered magne- 
tization E{n) is the expectation value ((— 1)'^'^"^), where 
s{n) is the number of spinons to the left of rung n. The 
average is taken over the ground state of nonrelativistic 
fermions occupying the first r levels with wavefunctions 
'>Pk{n) = •\/2/Lsin (jrkn/ L). The result for r spinons is 



E{n) = detS'(n). 
The r X r matrix S{n) has elements 



(21) 



sgn(2; — n) ip* {x)'ijjj (x) dx, (22) 



where "ipiix) are the spinon wavefunctions of the occupied 
states: i — 1 . . .r. See Appendix^for a derivation. The 
numerical data and the theoretical curve for = i/2 — 1 
(3 spinons) are shown in Fig.|2|[b). 

The agreement between the theoretical curve and the 
numerical data can be further improved by taking into 
account the spin AS^ ~ ^1/2 carried by the domain 
walls, 



_ 1 i-lTEin) 



2 ^ 

k=l 



(23) 



and by figuring in the reduction of staggered magneti- 
zation by quantum fiuctuations. The leading effect is 
a virtual exchange of a singlet and triplet on neighbor- 
ing rungs. This process increases the energy by J (the 
strength of triplet repulsion), has the matrix element 
5 J/2 ^ J, and thus can be treated as a small pertur- 
bation. To a leading order in SJ/J we find a simple re- 
duction of the envelope E{n) by the factor 1 — {SJ/J)"^. 
See Appendix El for details. 
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FIG. 10: (Color online) Friedel oscillations of the local mag- 
netization (S^) for three values of average magnetization m. 
The lines are best fits to the theoretical curves 1241 . The fits 
were performed in the range no < n < L+l — no with no — 5. 



Close to the plateau when the spinon gas is dilute, the 
deviation of the magnetization density from 1/2 is ex- 
pected to be proportional to \H — Hc2\^^'^, in complete 
analogy to the magnon case l|18|) . However, the magneti- 
zation curve (Fig. 0)) evidently remains linear almost all 
the way to Hc2, with only a hint of an upturn right next 
to the plateau. A possible reason for this behavior could 
be the narrowness of the critical region near Hc2 where 
the spinons can be treated as noninteracting fermions. 
One argument in favor of this interpretation is the rela- 
tive smallness of the square-root term whose amplitude 
is TT^^ jAS'^p/^ |2to*|^/^. In comparison to the magnons, 
the spinons have a reduced spin AS^ and a smaller mass 
TO* (by a factor of 4 to leading order in 6 J) . We also note 
that magnetization curves of an easy-axis XXZ chainii 
show a similar trend: the square-root term near the Neel- 
ordered state is relatively small. 



3. < m < 1; a Luttinger liquid 

The gapless phase in the field range Hci < H < Hc2 
is a Luttinger liquidi^ with continuously varying critical 
indices. The critical properties of the ladder are expected 
to be similar to those of the easy-axis XXZ chain in the 
gapless regime ffmin < |-Hxxz| < -ffmax in which the com- 
pressibility exponent K decreases monotonically between 
1 (dilute magnons) and 1/4 (dilute spinons). 

The compressibility exponent can be determined by 
examining the Friedel oscillations in local magnetization 
{S^) induced by the presence of the ends.^'^ In a ladder 
of length L, the leading behavior of magnetization away 
from the ends is 



const 



a cos (27rmn + (3) 



(24) 



[Lsin(7rn/L)]^ ' 
where to is the concentration of magnons and /? is a phase 
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FIG. 11: The compressibility exponent 7^ as a function of 
magnetization m in the gapless phase between fields H^i and 
Hc2- The error bars reflect the range of K values obtained 
for different choices of the short-range cutoff no. 

shift. Fits of the numerical data to this form are shown 
in Fig. ^1 for a ladder of length L — 160. The extracted 
exponent K is shown in Fig. ^2 a-s a function of the 
magnetization m. While there is a weak dependence on 
the short-range cutoff no, the trend is consistent with a 
monotonic decrease of K from 1 to 1/4. 

IV. INFLUENCE OF ANISOTROPY 

Experimentally accessible spin systems are almost in- 
evitably anisotropic. While an anisotropy may be small 
numerically, it may induce qualitative changes in the sys- 
tem's behavior because it lowers the symmetry. This is 
particularly important in the vicinity of a critical point 
or phase: the presence of a symmetry-breaking term 
can change the nature of a phase transition or even 
completely eliminate it. While the physical causes of 
anisotropics can vary, close to a critical point their in- 
fluence on the system can be expressed in terms of a 
few relevant physical variables. A well-known example 
is the breaking of the axial 0(2) symmetry of a spin 
chain under an external magnetic field by the anisotropic 
Dzyaloshinskii-Moriya (DM) interaction. The action of 
the DM term is equivalent to that of a weak staggered 
field transverse to the applied onei^^ 

A staggered field is uniquely defined in bipartite anti- 
ferromagnets. A frustrated antifcrromagnet can be par- 
titioned into two sublattices in more than one way. Ac- 
cordingly, several staggered fields can be introduced in 
such casesi^ Three staggered fields are potentially rele- 
vant to our ladder: 

L 

['^-o(-l)"5^ + ho^D-^ + h^^i-irD-^] (25) 

n=l 

(The fields are labeled by their Fourier indices.) The 
most relevant perturbation at a phase transition is that 



which couples directly to an order parameter. In our case 
( Ji < J2 < J), the "condensate" at Hd has the wavevec- 
tor (0, tt). The staggered field /iott coupled to it breaks 
down the 0(2) symmetry completely. In its presence, the 
quantum phase transition at Hd is expected to become 
a crossover. The other two staggered fields couple to the 
square of the order parameter, generating an easy-axis 
anisotropy along either x or y directions and thus lower- 
ing the symmetry from 0(2) down to Z2. In the absence 
of the more relevant staggered field h^^, the critical point 
at Hci will be preserved, but the universality class is ex- 
pected to change from commensurate- incommensurate to 
Ising. 

Our numerical work is focused on the influence of 
the staggered field /loir- (In what follows we drop the 
wavevector index.) We choose its amplitude to be pro- 
portional to the magnitude of the external field, as hap- 
pens with the Dzyaloshinskii-Moriya interactiontiSii^ h = 
cH. The anisotropy coefhcient c is varied between and 
0.1. Fig. 1121 shows the magnetization m{H) and energy 
gap A(i/) for c = 0.03 for several lengths of the ladder 
L. 

The introduction of the transverse field destroys the 
magnetization plateaux: the spin projection is no 
longer conserved. However, the fate of the quantum crit- 
ical points is different for the transitions out of the in- 
teger (i?ci) and fractional (iJc2 and Hcs) magnetization 
plateaux. A complete lack of finite-size effects near Hd 
is a good indication that the magnon condensation point 
Hd has become a crossover. In contrast, near the points 
of spinon condensation, Hc2 and -ffc3, the energy gap A 
is still sensitive to the system size, indicating that the 
critical points survive the introduction of anisotropy. 



A. Anisotropy and magnon condensation 

As noted above, the introduction of even a weak stag- 
gered field {h = 0.03H) completely suppresses finite-size 
effects in ladders of length L = 50 and more. This is 
fully consistent with the scaling theory of Bose conden- 
sation in the presence of a symmetry-breaking transverse 
fLeldi2^ A finite transverse field h — cH generates a fi- 
nite correlation length ^ oc A"-'^ cx c~^/^. Evidently, for 
c — 0.03 we have ^ ^ 50, so that ladders with L > 50 are 
already in the thermodynamic limit. 

The magnetization and energy gap for several values of 
anisotropy c are shown in Figs. IT^ a) and (b). The data 
are in agreement with the scaling theoryiSi As shown in 
Fig. ICT c). the energy gap obeys the scaling law 

A{H - Hd, c) = c^'^'mH - Hd)c-^h- (26) 

Exactly at Hd the energy gap is proportional to c^^^, 
while the magnetization m{Hd) oc c^^^. 
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FIG. 12: (Color online) Finite-size effects on magnetization 
(a) and energy gap (b) in the presence of a small anisotropy 
c = 0.03. The solid line shows m{H) for c = and L = 150. 
The dotted lines mark the critical fields Hci, H^i, and Hcz- 
L is the length of the ladders. 



B. Anisotropy and spinon condensation 

In contrast, for the same value of anisotropy, strong 
finite-size effects are observed near the spinon condensa- 
tion points Hc2 and i?c3 (Fig. I12|l . This indicates that 
the quantum phase transitions survive the introduction 
of the staggered transverse field /iott . 

The more robust nature of the spinon transitions can 
be traced to the spontaneous breaking of a discrete lat- 
tice symmetry at the fractional plateau: the two ground 
states (Fig. ^ break any symmetry transformation that 
exchanges even and odd lattice rungs. The addition of 
the staggered field /iott keeps some of these symmetries 
intact (e.g., the translation by one lattice spacing), so 
that the Hamiltonian remains more symmetric than the 
ground states. In other words, the Z2 translational or- 
der remains in what used to be the fractional plateau. 
Restoration of the Z2 symmetry at Hc2 (and Hcs) re- 
quires a phase transition, whether the transverse field 
/iott is present or not. 

The difference can also be understood by looking at 
the effect of the staggered field on the elementary excita- 
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FIG. 13: (Color online) Magnetization (a) and energy gap (b) 
near Hci for several values of the anisotropy parameter c (see 
text), (c) Scaling of the energy gap near Hci- The dotted line 
shows the location of the critical field. 



tions of the integer and fractional plateaux. By coupling 
to operators I?+ and Z?^, the staggered field adds or 
subtracts angular momentum 1. As a result, it creates 
or destroys single magnons (spin 1) at Hd and Hd but 
pairs of spinous (spin 1/2) at 11^2 and H^'^. A plausible 
effective Hamiltonian for the spinous with low momenta 
near the critical point in a weak transverse field would 
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FIG. 14: (Color online) Finite-size scaling of the energy gap 
near the critical field Hc2 in the presence of finite anisotropy 
c = 0.03. The dotted line marks the location of Hc2 in the 
absence of anisotropy. 

be 

Hxxz = '^[{p'^ /2m-n)alap+ivp{ala^_p-a-pap)], (27) 
p 

where p is the spinon momentum, ii oc H — Hc2 is the 
chemical potential, and v cx /iot is a pairing field. The 
energy gap behaves as follows: 

A=(l^l ^ ifM<-^'^ (28) 

^ymv (2fi + mv ) it ji > mv . 

The spinon condensation in 1+1 dimensions is thus 
similar to the commensurate-incommensurate transi- 
tions in two-dimensional statistical mechanics. In the 
absence of the transverse field /iott, it belongs to the 
metal-insulator (Pokrovsky-Talapov) universality class. 
Switching on the field converts the transition to the Ising 
universality class. 

We have verified that, for a fixed anisotropy c, the 
finite-size scaling of the energy gap is consistent with the 
Ising transition in 1-1-1 dimensions: 

A{SH,L)^L-'f{SHL), (29) 



FIG. 15: (Color online) Energy gap as a function of the field 
for several values of anisotropy c in a ladder of fixed length 
L = 50. The dotted lines mark the critical fields Hd, -ffc2, 
and Hc3- 

where SH = H — Hc^ic). Note that the critical field 
depends on the anisotropy parameter c. The scaling for 
c ~ 0.03 is shown in Fig. 1141 

As the anisotropy constant c increases, the critical 
fields Hc2 and Hcz shift towards each other, see Fig. ^| 
The two fields merge and the ordered phase disappears 
at a modest value of the anisotropy c « 0.06. Above this 
value of c, the ground state is non-degenerate and does 
not break the translational symmetry for any value of H . 

Ideally, we would have liked to verify the universal 
scaling at the quantum critical point H ~ Hc2, c = 0, 
1/L = 0, just like we did at the magnon condensa- 
tion point Hci- In this case, one expects to observe 
a crossover from the Ising critical behavior to that of 
the Pokrovsky-Talapov class. As can be inferred from 
Eq. 128|l . the energy gap is expected to cross over from 
\H — Hc2\ on the gapped side to\H — Hc2\^^^ in the for- 
merly gapless region. We have not been able to observe 
this crossover. This failure may be related to the nar- 
rowness of the region where the spinous can be treated 
as noninteracting fermions with a quadratic energy dis- 
persion. (See Sec. lIIIR2l ') To observe the crossover, one 
probably needs to work with a very small anisotropy c 
(we went down to 5 x 10^^) and rather long ladders to 
avoid finite-size effects. 



V. DISCUSSION 

We have presented a model of a one-dimensional quan- 
tum antiferromagnet in external magnetic field. The sys- 
tem exhibits both integer and fractional magnetization 
plateaux. The quantum critical points at the ends of 
an integer plateau are well-described by the Bose con- 
densation of magnons. At low densities, the condensing 
magnons behave as hard-core bosons or, alternatively. 
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can be described as weakly interacting fermions. The in- 
troduction of a weak anisotropy fully breaking the 0(2) 
symmetry of the model replaces the quantum phase tran- 
sition with a crossover. In contrast, a fractional magne- 
tization plateau breaks a Z2 (Ising-like) symmetry of the 
lattice and ends in a condensation of spinous - domain 
walls in the Zi order parameter. Magnetization patterns 
at low spinon densities are explained by modeling the 
spinons as free fermions. The introduction of a weak 
anisotropy preserves the Z-i symmetry of the model. As 
a result, the quantum phase transition at the end of a 
fractional plateau survives. 

This difference between integer and fractional plateaux 
is expected to show up in the properties of quan- 
tum antiferromagnets in higher dimensions as well. In 
that respect, SrCu2 (603)2, a physical realization of the 
Shastry-Sutherland model, is a prominent candidate, 
with plateaux at m = 1/8, m = 1/4 and m = 1/3 
that spontaneously break the translational symmetry of 
the crystaliSL2& It is well established by now that there 
are significant Dzyaloshinskii-Moriya interactions in that 
compound--- and that a gap persists in the region be- 
tween the m = and m = 1/8 plateaux, and closes 
(or has a deep minimum) before entering the m = 1/8 
plateau. This behavior is reminiscent of the gap behavior 
we found for the frustrated ladder. Given the pecularities 
of the triplet kinetic energy in the Shastry-Sutherland 
modelf^S*^ the extension to that case of the ideas devel- 
oped in the present paper is far from trivial. 
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APPENDIX A: THE SPIN OPERATORS ACTING 
ON A RUNG 



There are four states on each rung, which decompose 
into a singlet and a triplet under spin S = Si + S2. We 
denote the singlet as |s), and the Sz = 1,0,-1 compo- 
nents of the triplet as |-f), |0), and |— ) respectively. The 
spin difference operator D = Si — S2 acts as follows: 



D'-\s) = |0), £1^10) 
i?±|T)=±V2|,s), 



=±\/2|±), 



L»±|0) = 



(Al) 



APPENDIX B: SPINON ENERGY AT 0{5j'^) 

To compute the second-order correction to the XXZ 
Hamiltonian, one needs the action of Hi on states com- 
prised of singlets and 6*2 = 1 triplets. It follows from Eq. 
(jAlj) that D„ • D„+i acts non-trivially only on nearest- 
neighbor pairs of singlets, giving for one pair 



- 100)) 



(Bl) 



One can now do second-order perturbation theory to 
compute the correction to the XXZ Hamiltonian com- 
ing from annihilating and creating pairs of triplets out 
of the adjacent singlet states. Since these pairs are cre- 
ated between other magnons when the system is near 
the m = 1/2 plateau, the perturbation theory requires 
the diagonalization of "Hi on a length-4 chain of triplets. 
Calling the resultant states V'l a-nd their Tii eigenvalues 
Ai , one finds that the shift in energy coming from having 
singlets on adjacent sites is 



1(^.1^2 



C( J) 5 J2 



(B2) 



where 



C = 



1/20 1/5 (5-V2T)/20 (5 + V2T)/20 



-2 + J 



1+V2T 



J 



(B3) 

At J = 0.625 we obtain C = -2.026. The largest contri- 
bution comes from the last term where the excited state 
lies rather close to the spinon (an energy difference of 
0.255). 



APPENDIX C: MAGNON ENERGY AT 0{&.]'^ 



where D 



D"" ±iDy. 



In Sec. nil BTI we have determined the inverse magnon 
mass to leading order in 5J . Here we evaluate the next- 
order correction in the limit of well-separated energy 
scales Q. To this end, we consider the state of the lad- 
der with a single — ^\ triplet in the background of 
singlet states. For (5J = any such state is an eigen- 
state of the Hamiltonian Tio + Hi Q- The first-order 
correction in 5 J induces the hopping of the magnon to 
the next site with the amplitude (5 J/2, yielding an inverse 
mass 1/m* = \6J\. 

Second-order corrections to the magnon kinetic en- 
ergy include hopping through intermediate states with 
additional magnons. An example of such a process is 
shown in Fig. |31 Such a process involves the creation of 
a pair of magnons with total spin next to the original 
magnon and the subsequent destruction of a magnon pair 
by the perturbation term TLi ■ Since both of these events 
have the amplitude 0(5 J), we may treat the dynamics of 
the intermediate 3-magnon state using the zeroth-order 
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Hamiltonian TCq + 7^i. At that level, the 3-magnon com- 
posite is immobile. However, it has some internal dynam- 
ics: the spins of the individual magnons may fluctu- 
ate. We therefore first discuss the internal dynamics of 
the composite. 

A 3-magnon composite with the total spin = +1 
has 6 internal states: 



-00), |0+0), |00- 



In this basis. Ho — 3 — H, while 



Hi = J 



(CI) 
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(C2) 



The lowest energy composite state is close to that shown 
in Fig. I3b) and has the energy 3 — H — 3 J. The energy 
gap between the magnon and the 3-magnon composite 
is Ai = 2 — 3J. Our numerical studies were done on a 
system with J — 0.625, which yields a gap Ai = 0.125, 
comparable to the perturbation strength 6 J — 0.15. It is 
therefore not surprising that the leading-order result for 
the magnon mass was not reliable. 

To evaluate corrections at 0{6J^), we examine the 
hybridization of the magnon a,i with the symmetric 3- 
magnon composites (where 7 = 1 ... 6 is the internal 
index) : 



Wmagnon = ^((5 J/2) (aj,a„_i 4- H.c.) 



V C^, ^ C'y 



(C3) 



+ 5^ 5^ ^■y'^'yn^in 
n 7 — 1 
6 

+ ^^[4„(A^a„+i + p^Qn-l) + H-c], 
n 7—1 

where and pj are matrix elements involved in creating 
a 3-magnon composite centered on the left or on the right 
of the magnon: 

^ {s\cyn'H2al+i\s), /O7 = (s|c^„H2al_i|s), (C4) 

where \s) is the all-singlet vacuum state. Elimination of 
the hybridization term to 0{6J) via a unitary transfor- 
mation 



yn 



(C5) 



generates an additional magnon hopping term at 0{dJ^ 

6 



E^E«tia«+i+H.c. (C6) 

7—1 n 



This yields the inverse magnon mass 

p. 



7=1 ~' 

\5J\ + 2(\- ^ 



(C7) 



6 2- J 3(2 -3J) 



5J' 



The calculation is simplified because the one-magnon 
state only couples to three of the six composite states. 
For our choice of the coupling constants, the second-order 
correction to the inverse mass (0.575) is almost four times 
as large as the first-order term (0.15) thanks to the small 
energy gap between one- and 3-magnon states. 

The same unitary transformation allows us to deter- 
mine the shift of the magnon energy at this order: 



Aeia 



3(2 - 3.7) 



(C8) 



The second term comes from the effect of the magnon on 
the vacuum. The Hamiltonian term 7^2 creates virtual 
excitations in the form of two triplets on adjacent rungs. 
These virtual processes shift the vacuum energy by 



Ai.o = E^"'^'^^'' 



-3L5.P 



2 + 2J 4(2 -2 J) 



(C9) 



where \pi) is the intermediate state (jBl|l with two triplets 
next to each other. In the presence of a magnon the vac- 
uum fluctuations are suppressed in the immediate vicin- 
ity of the magnon. The factor L in Eq. (|C9I) is replaced 
with L — 4, increasing the energy of the 1-magnon state 
by the Casimir term in Eq. 



Casimir 



3(5JV2(1 - J). 



(CIO) 



The second-order correction to the magnon condensa- 
tion field is then 



AiJ, 



Aei. 



-9.5(5J^ 



-0.214. 



(Cll) 



Again, because of the small energy gap to excited states, 
the second- order correction to Hd exceeds the first-order 
one (—0.15) and does not improve the agreement with the 
numerics. 



APPENDIX D: ENVELOPE OF STAGGERED 
MAGNETIZATION 

In this Section we derive the expression for the enve- 
lope of the staggered magnetization (|21|l . We work in 
the continuum limit and treat spinous as noninteracting 
fermions in a box with < a; < L. Since each spinon 
serves as a domain wall, the envelope E(x) is found by 
averaging the operator 



E{x) = J]^ sg\\{xi - x) 

i=l 



(Dl) 
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over the (fully antisymmetric) wavefunction of r spinons 
1 

*(a;i...x^) = — ^e{fc}[|Vfe.(a;»), (D2) 

{k} i=l 

where = £ki...kr is the fully antisymmetric tensor 

with £l2...r = +1- 

The averaging yields 

^' {k} {1} »=1 

with the matrix Sij{x) defined in Eq. (|22l) . It can be 
simplified by noting that £{fc}e{;} = (—1)^, where P is 
the permutation that maps {fc} into {I}. Shifting from 
the sum over {k} and {1} — P({fc}) to a sum over {fc} 
and P allows the sum over {fc} to be performed trivially 
to obtain 

r 

E{x)^Y.^-\f\{S,pXx)=detS{x). (D4) 

P i=l 

APPENDIX E: REDUCTION OF STAGGERED 
MAGNETIZATION 

The domain wall magnetization envelope will be re- 
duced by the presence of quantum fluctuations. The 
leading order effect is the creation and subsequent annihi- 
lation of a pair of domain walls due to the interchange of 



a singlet and a triplet on neighboring sites. The domain 
wall Hamiltonian for the XXZ chain is^ 

xic]+iCj-i+Cj+iCj^i+ H.c), 

where the c operators obey fermionic commutation rela- 
tions. We shall ignore the quartic terms in this Hamilto- 
nian in the approximation of low spinon density. 

Since these domain walls represent boundaries between 
the Neel states with = (-l)"/2 and (-1)"+V2, the 
expectation value of the spin on rung n in a state {ip) is 

(Sn) = ^-^{mnm (E2) 

in addition to the spin carried by the spinons themselves. 
Here En = {—iy^"\ where s{n) is the number of spinons 
to the left of rung n. 

The spinon-number changing terms in the Hamiltonian 
can be eliminated to lowest order in e = 6J/ J by the 
unitary transformation 

j 

We can write the eigenstates of as — e^\^p), 
where jV') is an eigenstate of e~^He^ with definite 
spinon number. We now have that 



(-1)" 



(l-e2)(^|£;„|^)+0(l/L) + 0(e3). 



(E4) 



Here we have used the relation 

CmEn = sgn(m - n) EnCm (E5) 

and its adjoint. Note that the first order term vanishes 
because the operator [B, E^] includes a net change in the 
number of spinons. The 0(1/L) piece arises from terms 



in the commutator proportional to the spinon density 
operator. The second order term is independent of the 
number of spinons. Hence, the leading effect of the quan- 
tum fluctuations is to reduce the value of the staggered 
magnetization by a factor of 1 — {SJ/J)^. 
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